Gaussian coordinate systems for the Kerr metric 
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We present the whole class of Gaussian coordinate systems for the Kerr metric. This is achieved 
through the uses of the relationship between Gaussian observers and the relativistic Hamilton-Jacobi 
equation. We analyze the completeness of this coordinate system. In the appendix we present the 
equivalent JEK formulation of General Relativity - the so-called quasi-Maxwellian equations - which 
acquires a simpler form in the Gaussian coordinate system. We show how this set of equations can 
be used to obtain the internal metric of the Schwazschild solution, as a simple example. We suggest 
that this path can be followed to the search of the internal Kerr metric. 
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I. INTRODUCTION 



The recognition that natural processes cannot be in- 
fluenced by any choice of a representation of the events 
occurring in spacetime conducted the covariance princi- 
ple to be assumed as one of the fundamentals of modern 
physics. This idea was explicitly used to build a theory of 
gravity by Einstein, General Relativity, which gave a step 
forward by invoking that the MMG ( Manifold Mapping 
Group) should be taken as an invariance principle of the 
theory. In practical uses however, one is obliged to select 
a particular language by choosing a special coordinate 
system to describe a given phenomenon. Among all the 
possible choices one can make - most of them dictated by 
symmetry of the given problem - there is a very special 
one that can bring us more physical insight about the 
problem to be treated, that is the Gaussian coordinate 
system or Synchronic coordinate system. In fact, this co- 
ordinate system was suggested by C. Gauss in his works 
about curves and surfaces. In Gaussian coordinate sys- 
tem (GCS), a foliation of spacetime is made in such a way 
that one separates space and time, as in pre-relativistic 
theories. The time-like world-line of an observer that is 
orthogonal to the 3-d space (which is identified to the 
co-moving system) is such that the proper time of such 
observer coincides with the coordinate time. 

In the standard procedure made by the father founders 
of relativistic cosmology, the Gaussian coordinate sys- 
tem appears closely related to the Cosmological Princi- 
ple: the large scale structure of the Universe behaves 
as being homogeneous and isotropic. The solutions of 
Einstein equations which satisfies this postulate possess 
a complete Gaussian coordinate system, in addition the 
Gaussian surface is a Cauchy surface. 

On the other hand, there are cosmological solutions 
which do not satisfy this postulate. However these solu- 



tions have undesirable properties, as for instance, closed 
timelikc curves. 

A similar approach was done for Godel's metric by one 
of us [l| . It was shown that the Gaussian coordinate sys- 
tem is limited and cannot be extended beyond a certain 
region, the domain of which depends only on the vor- 
ticity present in this geometry. Such inaccessible region 
defines a frontier for time-like geodesies, which prohibits 
the extension of the GCS into the whole manifold. 

Mathematically, a Gaussian coordinate system is con- 
structed by the definition of a hypersurface S = S{x^), 
which satisfies 



as ds 

„Uiy dS dx^ _ n 



(1) 
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where are coordinates lying on S. The Gaussian co- 
ordinates are given by x^ = (S,x^). Eq. ([1} imposes 
= 1 and .gO* = 0. 

The main purpose of this work is to exhibit a Gaus- 
sian coordinate system for the Kerr metric. The method 
we use is provided by the relativistic Hamilton-Jacobi 
formalism of canonical transformations. The key idea 
is to identify the principal Hamilton function with the 
proper time of a test particle in this geometry. A im- 
mediate application of this method can be done for other 
solutions of Einstein equation as Schwarzschild and Kerr- 
Newman. Finally, in appendix]^ we list the coordinate 
systems encountered in literature for Kerr metric includ- 
ing our Gaussian coordinate system and in appendix [B] 
we regain from the Quasi-Maxwellian formalism (JEK 
equations together with the evolution equations of the 
kinematical quantities) the Schwarzschild solution in a 
Gaussian coordinate system. This formalism is an equiv- 
alent way to obtain the results of General Relativity. 
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II. CANONICAL TRANSFORMATIONS 
FORMALISM 

The action for a free test particle with mass m in the 
presence of a gravitational field can be written as 



ds. 



(2) 



Assuming that only one extreme of the path is fixed and 
the particle is moving on a geodesic line, we make the 
variational principle and discover how S depends on the 
coordinates x^. Such functional form is 



SS = —mUaSx", (3) 
where Ua is the 4-vclocity. Defining the 4-momentum 



Pa 



dS 



it satisfies 



PaP = m 



(4) 



(5) 



Replacing Eq. (|4]) into ([5]), we find the relativistic 
Hamilton- Jacobi equation for a test particle in a grav- 
itational field given by 



^x^^ dx^ 



(6) 



III. EQUIVALENCE BETWEEN A GAUSSIAN 
SYSTEM AND THE HAMILTON-JACOBI 
EQUATION 

By a canonical transformation of coordinates of a phys- 
ical system, so that {q^p) goes to (Q, P) via a generating 
function Fi{q,Q,t), we have the following constraint 



dFi = J2^p^dq, - P,dQ,) + {K- H)dt. (7) 

i 

In the covariant form, we write it as 



dF2 = -padq'' -Qf'dPf,. 
Then the function F2 must satisfy 



(10) 



Pa 



dF2 

dF2 



(11) 



Comparing with Eq. (U), we identify F2{q,P) as the ac- 
tion S{q,P) = S{q,p{q, P)) of a test particle in a gravi- 
tational field. Therefore, 



dS 
dS 

dp: 



(12a) 



(12b) 



Assuming = —S/m, from the 0-component of Eq. 
(|12b|) . wc get 



m 



dS 
dPo' 



and then, 



(13) 



(14) 



The Hamilton- Jacobi equation can be obtained making 
the new Hamiltonian K constant {H — > K = const) and, 
as we know that Pq oc K, it is possible to incorporate 
this constant into S and it follows that S = S{q'^'',Pi). 
Hereupon, rewriting Eq. (|6]) we obtain 



ago ago 

dx^^ dx" 



1. 



(15) 



We derive partially Eq. ([T5|) with respect to Pi and use 
Eq. (fT2b)) to get 



9' 



dx^ dx'^ 



0. 



(16) 



Note that the coordinate system together with Eqs. 
([T5|) and ([T6|) define a Gaussian coordinate system. 



dFi {q^,Q^) = -Pc^dq^ + P^C^Q^ (8) 

where p„ = {H,-p) and = {Kdt/dQ°,-P). Using 
another generating function F2{q^, P,y), such that 

P2 = Pi - P/.Q^, (9) 

wc have 



IV. APPLICATIONS 

A. Gaussian system for the Schwarzschild metric 

As a simple exercise, we exhibit a Gaussian coordinate 
system found from the definition ([T|) for internal and ex- 
ternal Schwarzschild solution. This case is particularly 
simple due to the symmetries of such metric. 
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1. The external case 

Considering Schwarzschild external geometry de- 
scribed in the usual coordinate system, for radial ob- 
servers with non null velocity at the infinity , the met- 
ric which is originally given by 

ds^ = (l - ^) dt" - (l - ^) dr^ - r^dn\ (17) 

where rn = 2M and M is the geometrical mass of the 
gravitational source, becomes 



transformation of (|T9|) if we choose — >■ 1 before the 
calculations, because we lose the degree of freedom nec- 
essary to build the other coordinates. In both cases we 
observe that the horizon "disappears", i.e., this coordi- 
nate system does not have any problem for r = r^j as 
in Eddington-Finkelstein or Kruskal-Szekeres coordinates 
system. However, in these cases, the natural "observers" 
arc null-type. On the other hand, in the Gaussian sys- 
tems (|19p and ((22)) the true natural observers ~ Sq 
are time- like in every point, making possible the descrip- 
tion of events of the spacetime from geodetic massive 
particles at rest. 



ds^ = dr^ - (a^ -1 + ^) dR^ - r(T, Rfdn^, (18) 
according to the following coordinate transformation 



T ^ at + Fe{r, a), 
u — dr 



(19) 



where a € 5R is a external parameter and r is interpreted 
as the proper time. Substituting this proposal in Eq. ([1]) 
one can see that F^{r,a) = dF^/dr must satisfy 



F^{r,a) 



(20) 



for the new coordinate system to be admissible. 

Analyzing the particular case in which the velocity is 
zero Q, we obtain the coordinate transformation inte- 
grating the following geodesic equations parameterized 
by the proper time t 



dH 
d\ 



A' dt dr 



(21) 



0, 



where A = [1 — rn /r) and wc obtain 



In 



y^r/rn + 1 

\\prh^- 1| 
o i2/3 



2 



(22) 



The new metric is obtained just choosing = 1 in Eq. 
p8)) . It means that we are describing particles whose me- 
chanical energy is equal to the rest energy. From this we 
conclude that > 1. Note that this coordinate trans- 
formation does not converge to the inverse coordinate 



The internal case 

To construct a Gaussian system for the interior solu- 
tion, we can consider a spherical shell filled with a perfect 
fluid of energy density p = const, and pressure p = co- 
moving to = (e''/^, 0, 0, 0). In this case the line element 
is given by 



where 



e--'d<^ - ( 1 - — j dr^-r^dn\ (23) 




(24) 



and tq is the radius of the "star" and rc = i/ p. We also 
assume that ro < Tc to avoid singulatities at the coor- 
dinate system. Similar to the coordinate transformation 
we have 



T = at + Fi{r, a). 



R = 



dr 



(25) 



where Fl{r,a) = dFi/dr must satisfy 



(26) 



Therefore, the new line element for Schwarzschild in- 
terior solution, written in Gaussian coordinates, is 



ds^ = dr^ - [a^ - e'')dR^ - r(r, R)^dn'^ 



(27) 



Of course, at r = ro Eqs. p7|) and ([25)) must be the 
same. Then we obtain that rn = r^/r'^. The neces- 
sary continuity condition of the metric on the star shell, 
in Schwarzschild coordinates, is naturally carried to the 
Gaussian coordinates. 
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Kerr Metric 



Substituting Eq. (|3T|) in Eq. ([T5|) . for m = 1, we obtain 



In 1963 Kerr found an exact solution jj] of Einstein 
equations, which describes the external space-time gener- 
ated by a "source" having geometrical mass M and angu- 
lar momentum a per unit of geometrical mass. In 1968, 
Carter developed a method to obtain all geodesies for 
this metric using the relativistic Hamilton- Jacobi equa- 
tion 0. In this paper, he considers the principal Hamil- 
ton function S* as a function of the proper time r of a 
test particle. Here we identify both. 

In the appendix we exhibit (see @), some known co- 
ordinate systems of the Kerr metric and their main char- 
acteristics. We include in this list the Gaussian system 
presented here. 

The line element for Kerr solution in Boycr-Lindquist 
coordinates is given by 



2Mr 



idr- 



p^de^ 



AM ra sin 

::2 — 



(28) 



(r2 



1^)8111^9 



2Mra'^ sin"' 



a? cos^ 9 and A 



2Mr. 



where = 

A simple way to calculate the geodesic equations pa- 
rameterized by the proper time for Kerr solution is to use 
the Euler-Langrange equations found in (3]. From first 
integrals of these equations we construct the tangent vec- 
tor field 



E 



2 



A (dWi 



p^A 



1 (dW2 

- sin^ C 
p^"Ash?y 



(32) 



■L'^ = 1, 



where E is interpreted as being the energy and L as be- 
ing the angular momentum for a test particle when we 
consider the asymptotically flat regime. We can rewrite 
this expressionin a more convenient form like that 



(dW2 



[EP ~ aLf 
Ea sin 9 — 



cos^ 9 = J, 



(33) 

where 7 is a constant (of separability). So, we have two 
equations, one for and other for W2 as 



Wi _ _^ v/(-Ef2-a£)2-A(r2+^ 



dr 



5 =±^/7-a2cos2 



(34) 



Ea sin 9 - J=^, 
sm t 



Note that dWi/dr = i^„(r) and dW2/d9 = G„(6'). 

The other coordinates, which we call here = R, 
= 9 and = $, are calculated from the spatial 
components of the second equation of ([T2)) by 



1/° -^[T.^E~2MraL], 



Fn{r) 



1 _ UIXLLA 
— 2 — 1\, 

P 

2 _ G„ {9) 



(29) 



1/3 [2Mr{aE - L csc^ 9) + p^L csc^ 9] , 

where = (r^ -I- a^)^ - a^Asin^ 9. The functions F„(r) 
and Gn{9) are given by 



_ J{EP^aLY-^{r^ + ^) 
^n[r) — ± ^ , 



Gn{9) =±Jj- a2 cos2 9 - (Ea sin 9 - -U 



(30) 



where = + a^ and E, L and 7 are constants of 
integration. 

Let us suppose that the new time coordinate is repre- 
sented by 5, a principal Hamilton function and written 
as follows 



R = 

e = 
$ = 



dS _,dWi _ dW2 

"5E~^ ~m' ~&E'' 



dS 
"5L 

dS 



, dWi dW2 /ocN 
UT' (^^^ 

91^1 _ dW2 



The remaining metric components can be given in the 
form 



511 

5l2 
513 



-EL 



2Mar sin^ 6* 



522 

523 
.933 



= -2[S(7-a2) + aL], 
= -{E^ + f"^ fAa^ 9) - 



2Mra'^ sin" 9 



-2 [L(r2 + 7) - {Ef'^ - aL)asm^ 9] 
-4:(r'^ + 7)(7 - cos^ 9). 
The determinant of this new metric is 



(36) 



S 



-Et + L(f> + Wi{r) + W2i9) = -mQ° 



(31) 



5 = detg^, = -4/(r)/i(0), 



(37) 
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where /(r) = [Ef"^ - aLf - A(r2 + 7) and h{e) = (7 - 
cos^ 9) sin^ 9 — (Ea sin^ 9 — L)^. If we guarantee that 
/, /i > in some region of Kerr spacetime, then the metric 
gfj,i, wiU be well-defined for these events. 

The usual coordinate systems for Kerr metric present 
problems in some regions, which are called horizons. In 
Boyer-Lindquist coordinates, for example, the mathe- 
matical expression for the Kerr horizons are 



{2E^ - l)sin2 



2EaL 



--^=-j=\^Er{E¥^ - ah) - rA - (r^ + 7)(r - Af)] , 

^ (40) 
where / and h were defined at Sec. [TVB] . So at the 
singularity r = and 9 = Tr/2 we obtain 



r± = A/ ± x/AP - a2 



(38) 



where ?'+(-) is called outer (inner) horizon. Differently, 
our Gaussian system presents a complete regularity at 
the metric components in the horizons. Note that there 
exists a divergence at the real singularity r = and 9 = 
tt/2. 



lim 



[Q - a\E^ ^ 1) + L""] 
cos 9^/Q 



lim 



7Af 



^0 ar' 



(41) 

where Q = 7 — {Ea — L)^ . We conclude that the diver- 
gence — ^ —00 at the singularity is guaranteed just for 
the congruence which have Q = 0. All the other congru- 
ences cannot reach the singularity. 

If we pick out a given black hole with A/ = 2 and = 



We choose a = 5^ in this coordinate system, which 1 (fixed that > these numbers are completely 



is obviously geodetic. Taking the inverse coordinate 
transformation, we write the observers field in Kerr co- 
ordinates in such way that 



arbitrary and they do not interfere at the results), we can 
make a more detailed analysis controlling the boundaries 
of the functions /(r) and h{9). Substituting x = sir? 9 
in h{9), we get 



y° = ^P^-B - 2MraL\, 



T/2 



A, 



(39) 



= ^ [2Mr{aE - L csc^ 9) -t- p'^L csc^ 9] . 

If we compare Eq. ([55)1 with Eq. ([^ obtained 
from the Euler-Langrange equations we conclude 
that they are the same. Therefore, 5q corresponds to 
all tangent vectors of timelike geodesies for Kerr solution. 



V. COMPLETENESS OF THE GAUSSIAN 
COORDINATE SYSTEM 

At present there is not a theorem that specifies the 
necessary and sufficient conditions of existence of a com- 
plete Gaussian coordinate system for an arbitrary metric 
or whether the hyper-surface defined by it is a Cauchy 
surface Q. 

In our case, it is obvious the impossibility to identify 
S with a Cauchy surface due to the presence of closed 
time-like curves. What about the completeness of such 
Gaussian system? To answer this question one can ana- 
lyze eventual divergences in the expansion factor i^. This 
will allow us to recognize the regions where a given con- 
gruence of observers remains well-defined or not. A given 
congruence of observers is nothing but a given choice of 
the parameters i?, L and 7. The explicit expression of 
the expansion is 



h{x) = (1 - E^)x^ + (7 - 1 + 2EL)x - 



(42) 



By the roots of this polynomial function, we can en- 
counter intervals of the parameters range such that we 
can cover all values of 9. Thus, we can choose L = and 
E^ > 1, and consequently we obtain E^ < 7. By the 
other hand the /(r) function becomes 

/(r) = {E^-l)r^ + 4:r^ + {2E^-l-jy + 4:-/r+{E^--/). 

(43) 

From the zero-order term in r, we see that E^ > 7 IS a 
necessary condition for f{r) to be greater than for all 
values of r. Therefore, we conclude that is impossible 
to cover all manifold events with only one congruence of 
observers. The case Q = provides the only congruence 
which can cover all values of 9 and all positive values of 
r. 



VI. GENERALIZATION (KERR-NEWMAN 
METRIC) 

In 1965 E. Newman et al. Q found a generalization of 
Kerr solution, which describes a black hole with geomet- 
rical mass Af, angular momentum a per unit of geomet- 
rical mass and charge q. 

The line element for Kcrr-Newman solution is 



= 1 _ 



2(2Afr-g2)asin^i 



(44) 



(r^ + cr) sm 9 + ^ ^ 

P 
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where = r"^ + cos^ Q and A = + + — 2Mr. 

If we make all steps used in the previous section, at the 
end we find the covariant metric components for Kerr- 
Newman solution described by a Gaussian coordinate 
system, which can be given as follows 



511 



1 



2Mr- (?2 



hL , 



.912 = 

313 = -2[S(7-a2) + aL], 



522 = -(L2+f2gij^2^^_ 



(45) 



{2Mr-q^)a^ sin* 6* 
7 



323 = -2 + 7) - {Ef'^ - aL)asm^ 



933 



-4(r^ + 7)(7 — cos^ 



Now = Sq corresponds to all tangent vectors of 
timelike geodesies for Kerr-Ncwman solution. 

If we take appropriate limits wc get Gaussian coor- 
dinate systems for other metrics: in the case a = 
and q we obtain a Gaussian system for Reissner- 
Nordstrm solution, if a 7^ and q = we construct the 
previous Gaussian system for the Kerr metric, and fi- 
nally, if a = g = we obtain the Gaussian system for 
Schwarzschild case. 



VII. CONCLUSIONS 

Following a systematic algorithm, we have built a set 
of Gaussian coordinate systems for the Kerr metric and 
generalizations. We analyzed the completeness of the 
Gaussian system which is intrinsically related to the con- 
served quantities associated to a test observer immersed 
in this geometry. Another important feature is the fact 
that in this coordinate system the K-metric is not static, 
i.e., a Gaussian observer measuring the geometric prop- 
erties in its neighborhoods, concludes that such metric 
is not static according to its own proper time. In the 
appendix we show how it is possible to use such GCS in 
order to undertake the search of the corresponding inter- 
nal solution. This work is under analysis. 



transformations between them. Here there is a list of 
them, with some properties, and the inclusion of our 
Gaussian system. 

• Kerr's original coordinates: 

In this coordinate system {u,r,0,(j)), the line ele- 
ment is 



= 1 - 



2Mr \ 



"2— y du^ - p^dd'^ - 2as\T? 9drd(j) 



^2dudr ~ ^MiapLAdudc^^ 



(r2 + a") sin^ 9 + 2Af ra^ sin^ 
P 



(Al) 



where = + a} cos^ B. 

The most important features in this coordinate sys- 
tem are: 

— The appearance of a actual singularity in gum 
for r = and 6* = 7r/2; 

— Setting a — 0, the line element re- 
duces to Schwarzschild geometry in "advanced 
Eddington-Finkelstein coordinates" ; 

— In terms of M, the line element can be 



put into Kerr-Schild form by g 



9o 



f{M,a,x°')m'' , where g'^" is Minkowski 
spacetime and a geodetic null vector. 

• Kerr-Schild "Cartesian" coordinates: 

By use of (a;°, x, j/, z) coordinates, the ds^ for Kerr 
metric becomes 



ds^ = {dx^Y - dx'^ - dif - dz^+ 



dx^ + ^2 j^^2 {xdx + ydy)+ (^2) 
2 



'i'^^i iydx-xdy) + ^dz 



The main characteristics of this coordinate system 
are: 
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Appendix A: Some Special coordinate systems for 
the Kerr metric 



As we said, in [6| we see many coordinate systems for 
Kerr metric encountered in literature and the coordinate 



— For M — > 0, it is Minkowski spacetime in 
Cartesian coordinates; 

— For a — >■ 0, it is Schwarzschild solution in 
Cartesian coordinates; 

— The full metric {M, a 7^ 0) is obviously the 
Kerr-Schild form again. 

• Boyer-Lindquist coordinates: 

The most useful coordinate system for Kerr metric 
is Boyer-Lindquist coordinates, as follows 
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ds^^il- j dp -^dr^- p^dO^ 
^2 , 2Mar -2 , 



(A3) 



where A = + — 2Mr. 

The noticeable properties are that: 

— It minimizes the number of off-diagonal com- 
ponents of the metric; 

— The asymptotic behaviour of this coordinates 
permits to conclude that M is indeed the mass 
and J = Ma is the angular momentum; 

— For a —> 0, it is Schwarzschild solution in stan- 
dard coordinate system; 

— For M — > 0, it is Minkowski line element in 
oblate spheroidal coordinates; 

— It is a maximal extension of the Kerr manifold. 

• Rational polynomial coordinates: 

If we make a coordinate transformation x = cos Q 
from Boyer-Lindiquist coordinates, we have the fol- 
lowing new version for Kerr spacetime 



— For a 0, it is Schwarzschild geometry in 
Painlev-Gullstrand form; 

— For M — > 0, we obtain Minkowski spacetime 
in oblate spheroidal coordinates; 

— In Doran coordinates, the contravariant met- 
ric component g*"^ is equal to 1; 

— According to ADM formalism, Doran coor- 
dinates slice the Kerr metric such that the 
"lapse" function is everywhere unity. 

• Gaussian coordinates: 

Constructed from the relativistic Hamilton- Jacobi 
equation. The line element for the Kerr Metric in 
a Gaussian coordinate system is 



2Mr \ 



dB? 



i2 B + 2MraW0 \ ^qs^ 



^4(7.2 + - cos^ 6')(i$2 + 
EL + 2Marsin2g ^ dMe+ 

-2 [L(r2 + 7) - [Ef"^ - aL)a sin^ 6] dOd^, 



ds^ 



1 



2Mr 



dP- 



r-t^dr^ 



where 



1- X T + a. X 



(A4) 



-(1 - x ) + H ■> , '2 i 



dcl,^ 



These coordinates introduce the following qualities: 

— All metric components are rational polinomial 
of the coordinates; 

— The non-appearance of trigonometric func- 
tions do the computational calculations faster 

• Doran coordinates: 

Introduced by C. Doran in 2000, here we obtain 
another coordinate system for Kerr metric given 
by 



ds^ - dp - p^de^ - (r2 + a") sin2 Odc^^ 



(dt — a sin 



(A5) 



The useful features of Doran coordinates are: 



Besides the usual characteristics of the Gaussian 
coordinates, in this case we have that 

— For massive test particles, the geodesic equa- 
tions parameterized by the proper time are im- 
mediately integrated; 

— Differently from all other cases, the metric is 
non-static but it is stationary; 

— This metric depends on parameters of the ob- 
servers field comoving to the reference frame; 



Appendix B: Schwarzschild solution from 
quasi-Maxwellian equations written in Gaussian 
coordinates 

As it was showed in [lOi, the JEK (Jordan-Ehlers- 
Kundt) equations are equivalent to general relativity. Be- 
sides, they become particularly simple when expressed in 
a Gaussian coordinate system. Therefore, in principle, 
we could use this simplicity to search of an internal so- 
lution for the Kerr metric. As an example, we will apply 
this method to obtain the stellar Schwarzschild solution. 

There are many references treating the formal deduc- 
tion of these equations and their properties as for in- 
stance in [ll|. These equations (JEK) can be obtained 
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from Bianchi's identities, together with Einstein equa- 
tion, that is 



(Bl) 



From this we obtain the corresponding independent pro- 
jections of the divergence of Weyl tensor 



(B2) 



Besides, we also have the conservation law of the 
energy-momentum tensor 



(B3) 



which can be projected parallel to an observer field 
or perpendicular to it, as it follows 



(B4) 



where h^'" = g^"" -V^'V. 

From the Riemann tensor, we can find the evolution 
equations for the kinematical quantities expansion -d, 
shear a^i^ and vorticity given by 



1? 



1 2 

1 9 



together with the constraints 



(B5) 



+ 2a;"aa = 0, 

-^/l(r^/iA)"?7e'''^''K(a'Q/3 + ^^a/3);7 + a(TWA) = -H"rA- 

(B6) 

If we assume that Einstein equation is only valid in 
a Cauchy surface, the set of equations (|B2p - (jB6p . so- 
called the quasi- Maxwellian equations, propagates it to 
the whole spacetime. 

Let us consider a diagonal metric, similar to 
Schwazschild one described in Gaussian coordinates, as 
follows 



ds-" 



dT'^ - B{T, R)dR^ - r2(T, R)d^f 



(B7) 



and an observer field = 5^. The expansion d for this 
vector is given by 



(B8) 



where Y{T, R) + dY/dT. After that, we calculate the 
shear tensor a'^i, and the eletric part of Weyl tensor 
{Ef^i, = —Wa^/SuV^V^) and write them in matricial form 





(B9) 



1 



[E^,]^g{T,R) 



where 



f(T,R)^ 



0- 



l B 2f 



1 



(BIO) 



3 \ B 



(Bll) 



and 

9{T,R) = 



1 



12r^B- 



+2rBfB + 2rr'B' 
+4Br'^), 



■{-2r'^BB + r^B'^ - ArBr" 
ArB'^r - AB'^ - 4B'^P + 



(B12) 



where Y'{T, R) = dY/dR. 

Observe that both cr^j^ and -E^j/ are proportional. All 
other quantities like the magnetic part of Weyl tensor 
{Hap —*Wafj.p,yV^V'^), vorticity ujap and acceleration 
(a^ = V^^-^V) are identically zero, due to properties of 
the observers congruence chosen 

Let us assume that 
co-moving to an arbitrary fluid, which can be expressed 

by 



Sq in Gaussian coordinates is 



Tf,u = {p + p)Vf_y^ - pgt^i, 



(B13) 



where p is the energy density, p is the isotropic pressure, 
is the heat flux and tt^^ is the anisotropic pressure. 
In the case of Schwarzschild stellar solution, it is as- 
sumed a perfect fluid inside a spherical shell and an 
accelerated observer [u^ = -y/ffoo "J^) comoving to this 
fluid. With respect to the Gaussian observers such fluid 
presents a heat flux = (0,q^,0, 0) and an anisotropic 
pressure 



9 



/ 1 



n{T,R) 









1 

"2 



(B14) 



With these considerations, the set of equations (jB2[) 
takes the form 



<7 + 3fg = \fTi-\{p + p)f+\TT + ldn+ (B15) 



1 B' 



The conservation law (|B4p can be written exphcitly as 



p+(p + p)i5-|/7r+(gi)'+(l^ + 2^) =0, 
tt' + 3^7r + 91 + - p' = 0. 

(B16) 

Now, the evolution of the kineniatical quantities is the 
following 



where h{R), y{R) and b{R) are arbitrary functions and 
fc is a constant. If we assume as initial condition sur- 
face r{T,R) = const — > 00, then we will obtain the 
Schwarzschild external solution. 

The Schwarzschild internal solution can be obtained if 
we consider the energy-momentum tensor associated to 
the Gaussian observer Sq , written in terms of the quanti- 
ties associated to the observer (energy density p and 
pressure p) presented in Sec. |I V A 2j . as follows 



Pg = (p + p)a'^e " -p 
PG = ~\[{p + p){l-a^e-'')-ip] 
g^ = (p+p)ae-''(0, 1,0,0) 
7r= hl-a^e-") 



(B20) 



Finally, the only constraint equation is 

/' + 3-/ - = 0. (B18) 
r 6 

Hereupon the set of equations (jB15P - (jB18[) corresponds 
to a problem of initial conditions, which shall give origin 
to Schwarzschild solution. As it is not our aim, we will 
not solve these equations step by step. However, we will 
indicate how to proceed. 

First, we analyze Einstein equations for T^^, = and 
we obviously obtain B{T, R) and r{T, R) like those given 
in Eq. p8|) . identifying T = t. On the other hand, from 
the quasi-Maxwellian equations we get 

f = ^y{R) + ^, (B19) 
/ = b{R)^Jh{R) + ^ 
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